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RESULTS OF INCLUDING A BOUNDARY SHOCK WAVE IN THE 


CALCULATION OF THE FLIGHT PARAMETERS OF 
A LARGE HIGH-SPEED METEOR 
By E. Dale Martin 
Ames Research Center 


SUMMARY 


The results of including a boundary shock wave in the calculation of the 
flight parameters of a large high-speed meteor with rapid vapor ablation are 
investigated. The boundary shock wave, a thin viscous region in a gas flowing 
rapidly (at high Reynolds number) out of a surface, is characterized by heat 
conduction and viscous -compressive stresses such as occur in a shock wave. 

The viscous effects represent the translational nonequilibrium induced in the 
extreme cases of very rapid vaporization produced by absorption of intense 
heat radiation at the surface. 


Calculations for the particular case of the Ondhejov meteor Pribram use a 
data curve fitted graphically to the meteor-tracking data* These calculations, 
based on the approximate tracking data and on gross approximate equations 
describing the vaporization ^process, indicate that the proper range of condi- 
tions was present on the Pribram meteor, for a significant portion of its 
trajectory, for a very thin, strong boundary shock to occur at the meteor 
surface in the vapor, according to the boundary- shock-' wave theory. The abla- 
tion Reynolds number was high, and significant ratios of density and pressure 
across the boundary shock wave result from the calculations • The meteor 
entered the atmosphere at about 20*9 km/sec* At an altitude of 55 km, where 
the meteor radius was indicated in the idealized calculation to be l8 .6 cm, 
the order of magnitude of the calculated boundary-shock thickness was in the 
order of l/lOOO times the entire thickness of the vapor layer, or about 
7X10" 4 cm. Although the usual convective heat transfer from the hot air layer 
was completely blocked by the efflux of vapor, an estimate according to this 
theory indicates high heat conduction at the surface in the boundary shock 
(at least 20 percent as large as the radiative heating). 

Accounting for the effects of the boundary shock wave substantially 
decreases the pressure drag, but the viscous-compressive drag is high, so the 
effective drag coefficient (at least in this simplified analysis) is virtually 
unchanged by the presence of the boundary shock. Although heat conduction in 
the boundary shock is high, it is balanced by the work of blowing off the 
vapor against the resistive force due to the viscous stress, and the overall 
heating of the body is little affected. 

Detailed calculations of the flow about a large high-speed meteor would 
be significantly influenced by the boundary-shock effects, but overall gross 



values of the heating and motion parameters calculated by methods used 
recently by Hi J. Allen and N. A. James are valid because of their use of 
appropriate assumed values of certain parameters. 


INTRODUCTION 


An investigation is undertaken to determine what possible effects a 
boundary shock wave may have in a simplified calculation of the flight param- 
eters of a large high-speed meteor and to determine qualitatively the signifi- 
cance of including consideration of the boundary shock wave in this 
calculation. A recent theoretical study (ref. l) investigated the possible 
occurrence of the boundary shock wave, a thin region where viscous-flow 
effects predominate in the rapidly ablating vapor very near the molten surface 
of a body, such as a meteor under certain conditions, moving at a very high 
speed through the atmosphere. The viscous effects in the vapor are effects 
of the translational nonequilibrium induced by the very high rate of vaporiza- 
tion with its attendant high rate of heat conduction into the body. 1 In the 
thin viscous region at the body surface in the vapor, viscous -compressive 
stresses (as in a shock wave), accompanied by heat conduction, would be 
present, in contrast to the viscous shearing stresses present in boundary- 
layer-type flows • However, the layer is much like a boundary layer in that it 
is adjacent to a wall and has inviscid flow on only one side. The thickness 
of the viscous region becomes vanishingly small as the Reynolds number of the 
efflux becomes large, as is true for either a shock wave or a boundary layer 
when the appropriate Reynolds number is large. 

The rapid change of the flow variables across a boundary shock wave could 
significantly affect the flow pattern of the vapor in front of the meteor. 

One would then be interested in whether the viscous effects, if they occur, 
might also influence significantly the motion and heating of the meteor. The 
purpose here, then, is to investigate a known physical situation having condi- 
tions under which the boundary shock wave would be possible (according to the 
theory of ref. l) in order to determine the effects. 

The importance to aerodynamic ists , who^are concerned with hypervelocity 
continuum flow, of acquiring and being able to understand and interpret the 
data from flight observations of meteors has been discussed recently by 
Allen (ref. 2 ) and Allen and James (ref. 3). Before knowledge of meteor 
phenomena can be applied to the technology of high-speed missiles and space 
vehicles, aerodynamic ists must have sufficient understanding to predict at 
least approximately the motion and heating of meteoric bodies with known 
characteristics. Recent analytical studies (see ref. 3 ) have been very 
successful in this respect. 

In the analysis presented here, attention will be confined to stone 
meteors with vapor ablation rapid enough that the vapor flow divides into 


lA more detailed description of boundary shock waves is contained in the 
section immediately following the Introduction. 
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inviscid and thin viscous regions (see Analysis)* Of particular interest is 
the Ondrejov meteor Pribram, a stone meteorite that fell in Czechoslovakia in 
1959* for which tracking data were obtained by the Ondrejov Observatory 
(ref. k) and for which the mass density is known from recovered fragments. 

The Pribram meteor undoubtedly had a high rate of vapor ablation over most of 
the portion of the trajectory for which the data were taken, so it is a good 
example for the study of the boundary-shock effects. 

Many idealizations are made in the analysis because of unknown physical 
properties, but the order- of -magnitude results obtained are expected to be 
qualitatively useful. Precise calculation of the motion and heating of a 
large high-speed meteor to determine the effects of a boundary shock wave 
would require detailed analysis of the flow field surrounding the meteor. 
However, it was deemed appropriate to approach the calculation here in much 
the same manner as in reference 3* Rather than a detailed flow-field analysis, 
then, this study uses the meteor-tracking data for velocity and acceleration, 
which are carefully fitted graphically to what is believed to be the most 
reasonable form for the data curve. For the calculation, a set of equations 
to determine the appropriate meteor parameters is formulated, which includes 
use of: (l) the previously derived results of reference 1 for the equations 

relating conditions across the boundary shock wave, (2) conservation equations 
and other relations incorporating certain appropriate assumptions about the 
flow around the meteor, and (3) points from the curve fitted to the meteor- 
tracking data as input. The approximate analysis assumes one -dimensional 
flow of the ablating vapor and constant average values of flow variables over 
the meteor face and employs a shape factor to account for either a sphere or a 
flat-face body. For gross results, such an approach is customary (ref. 3) 
except that the shape factor was not needed before (see below). More complete 
conservation equations than have been used previously are derived. These 
equations and the derivation contain useful information not obtainable from 
less complete statements of the conservation principles. 


THE NATURE OF, AND CONDITIONS FOR OCCURRENCE OF, A BOUNDARY 
SHOCK WAVE IN RAPID VAPORIZATION 


A boundary shock wave is most easily understood as a region of 
translational nonequilibrium in a molecular flow at a very high rate out of a 
surface. (Adequate understanding of it therefore requires some background in 
the basic concepts of modern kinetic theory, for which the reader is referred, 
e.g., to refs. 5"7 • ) 

Translational nonequilbrium in a molecular flow is manifested 
macroscopic ally by the "transport processes" of viscosity and heat conduction . 
Translational nonequilibrium is also induced by application of significant 
viscous stress or heat conduction in, or at the boundaries of, any region of 
gas flow. (The familiar fluid-dynamic boundary layer is a region of transla- 
tional nonequilibrium induced by shear at a wall; the familiar gas -dynamic 
shock wave can be thought of as a region of translational nonequilibrium 
induced by the diffusion or mixing, and subsequent rapid accommodation, of the 
molecules of one equilibrium average state into a different equilibrium 
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average state, which necessarily entails heat conduction.) Sufficiently far 
from any source of translational nonequilibrium, the gas flow relaxes to a 
local-equilibrium state (or decays to an "inviscid f low" ) • "Sufficiently far 
from any source . . ."is equivalent to saying "for a sufficiently large 
Reynolds number based on distance from the source" or to "for a sufficiently 
large distance for diffusion of the transport phenomena" (cf. ref. 8, 
especially pp. 21 - 23 ) • 

During vaporization at a very high rate, produced by absorption of 
intense heat radiation and subsequent heat conduction back toward the liquid 
in the vaporization region, it is expected that as the vaporization proceeds 
to a suff iciently high rate, the associated heat conduction will become strong 
enough to induce significant translational nonequilibrium in the vapor as it 
leaves the molten surface. The translational nonequilibrium is induced 
because of the lag in transferring energy from the molecular translational 
degree of freedom normal to the surface to the lateral degrees of freedom of 
the molecules. The lag occurs because of the spaces between the molecules as 
they are separating in the phase-change region. From kinetic theory, this 
lag, or difference, in the distribution of molecular energy in one direction 
from that in another direction is understood as the essence of translational 
nonequilibrium, including viscosity and heat conduction. This translational 
nonequilibrium (viscosity and heat conduction) then requires a certain number 
of collisions, and hence a certain time and distance of flow from the wall, to 
"equilibrate,” or to relax to a condition of flow in local translational 
equilibrium (inviscid flow). The relaxation distance is the thickness of the 
boundary shock wave . 

From the Chapman- Enskog procedure of kinetic theory, it is found that a 
gas flow in local translational equilibrium is governed by the Euler equations 
of inviscid flow and that regions of flow not too far from local translational 
equilibrium are adequately described by the Wavier- Stokes equations of viscous 
flow. It is known (e.g., ref. 9 ) that the structure of a sufficiently weak 
shock wave (and of a downstream portion of a strong shock wave) is accurately 
described by a solution of the Eavier-Stokes equations. 

In a boundary shock wave the viscous stress is compressive (rather than 
shearing) as in a normal shock wave; and also, as in a shock wave, the viscous 
dissipation is accompanied by heat conduction back through the gas flow. The 
same equations govern the flow through a plane boundary shock (or through a 
very thin curved boundary shock) as govern the flow through a shock wave. 
Reference 1 used the Eavier-Stokes formulation to obtain the conditions across 
a boundary shock, which are represented by equations given below in the sub- 
section "Equations for Ablation Process and Boundary Shock Wave." In a pri- 
vate communication. Prof. H. Schlichting has referred to the solution for the 
boundary shock wave given in reference 1 as "a . . . generalization of the 
known solution for the normal shock wave." It is a generalization in the 
sense that a boundary at upstream infinity can be brought into the region of 
the rapid transition; or alternatively, the limiting, or degenerate, case of 
a boundary shock wave as the boundary moves upstream to infinity relative to 
the rapid -transit ion region, or equivalently as the boundary shock "detaches" 
from the wall (as the heat conduction vanishes in a supersonic efflux), is a 
simple shock wave. From another point of view, a wall can simply be inserted 
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at any point in a simple shock-structure solution; the values of the various 
flow quantities at that point are then the boundary conditions , and the down- 
stream portion of the original shock structure then represents the structure 
of a possible boundary shock wave. 

In reference 10, the necessary and sufficient conditions for a boundary 
shock wave to occur were shown to be 

(a) large Reynolds number based on conditions at the boundary of 
the vapor , and 

(b) significant heat conduction or viscous stress (significant 
translational nonequilibrium) at the boundary. 

More thorough discussion of boundary shock waves and the development of the 
theory is contained in reference 10, including description of the asymptotic 
treatment of the translational nonequilibrium flow, with emphasis on the 
boundary conditions, and including discussion of the appropriateness of the 
equations used in the approximate treatment of the phase change. 


ANALYSIS 


General 

A large meteoric body of the type that would experience a high rate of 
vapor ablation while moving at very high speed through the earth 1 s atmosphere 
is shown in sketch (a). 2 If the observer were moving with the meteor, one 
would expect the instantaneous flow pattern to be qualitatively as shown in 
the sketch. (The curved lines, except the one indicated as the shock wave, 
are streamlines . ) The very high rate of vapor ablation is to be expected for 
a large high-speed stone meteor such as the Ondrejov meteor Pribram (see 
ref. 3). 


2 The meteoric body sketched is nearly spherical. Very large, very high- 
speed meteors would be expected to have a nearly spherical shape because of 
their high rates of ablation and their general rotational motion. Virtually 
no spherical meteorites are found, however, for the following reasons: 
Meteoric bodies that are either too small or too slow will not have experi- 
enced the extreme ablation and hence will, in general, be irregular in shape 
(see ref. ll). Nearly all (if not all) the fireballs that reach the earth 
without vaporizing completely will have fragmented at some point of the 
trajectory (ref. 11 ) and hence will also be irregular. 
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Shock wave 



Sketch (a).- Flow configurati 


Although a detailed flow-field analysis is not made here, it is important 
to understand the qualitative structure of the flow. The following character- 
istics of the structure should he noted (refer to both sketches (a) and (b)): 

A stagnation point exists at point 0 , on the interface between the hot com- 
pressed air and the vapor ablating from the meteor face. For large Reynolds 
number p2U2d 0 /p2 (notation is defined in appendix A), the flow of hot com- 
pressed air around the meteor is essentially inviscid except for the very thin 
shock wave and a thin layer at the interface. For large Reynolds number, 
Pb u bdi/ Mb*> meteoric vapor is also essentially inviscid except for 

the thin viscous region at the interface and the thin boundary -shock region 
(see ref. l). The three viscous regions - the shock wave, the interfacial 
layer, and the boundary shock wave - may be treated as surfaces of discontinu- 
ity of certain variables in the inviscid flow. It can be shown that only the 
pressure is required to be continuous across the interface (in the inviscid 
limit). Hence, the pressure p Q is the stagnation pressure on the axial 
streamline between point 2 and point e in sketch (b). 

The high temperature air in the outer shock layer radiates heat to the 
meteor surface at a high rate. The large energy flux to the surface produces 
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Sketch (h).- Viscous and inviscid regions in flow near body axis. 

the high rate of vapor ablation . 3 The hot air and vapor are arbitrarily 
assumed to he transparent to radiation, and the radiation is assumed to he 

3 For small or relatively slow meteors, heating hy friction with the air 
(convective heating) produces ablation of the meteor surface. In the case of 
the fireball, at least in the major portion of the trajectory where the heat- 
ing is most intense, the meteor is being heated not by friction with the air 
but by heat radiation from the highly compressed air behind the shock wave. 
Heating by friction, or convective heat transfer to the body, is almost non- 
existent in this case, as will be seen later. However, an additional source 
of heat conduction is considered here: the viscous dissipation associated 

with the translational nonequilibrium in the thin region adjacent to the 
liquid layer. 
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absorbed -within a very short distance inside the liquid surface. Radiation 
absorption by the vapor may be important. However, 

(a) The radiation properties of the vapor are unknown; therefore, some 
assumption must be made in order to study the qualitative features of the 
flow; 


(b) Neglect of the absorption by the vapor simply changes the magnitude 
of the radiation absorbed at the liquid surface, which should not alter the 
qualitative results unless the amount of radiation received by the body is 
changed by orders of magnitude; 

(c) The vapor is at a much lower temperature than the hot air in the 
outer layer, which tends to decrease the importance of absorption by the vapor. 

Since the meteoric -stone material has low thermal conductivity, only the 
surface is at the vaporization temperature; it is assumed that vapor forms 
only at the surface, and little or no bubble formation is expected to accom- 
pany the "boiling off” of the vapor. In addition, the high viscosity of the 
molten meteoric stone prevents the spraying of liquid drops along with the 
efflux of vapor (ref. 12). 

Allen (ref. 3) has pointed out that meteors of the type considered here 
tend to be spherical, but may be somewhat flattened by ablation due to large 
rates of radiative heating (see sketch (a)). Moreover, the present desire is 
to obtain only gross effects; we therefore take advantage of the great simpli- 
fication provided in the analysis by roughly approximating the flow as being 
one -dimensional near the meteor surface. That is, we take constant values for 
the variables over the meteor face and consider the area of the meteor face to 
be equal to the cross-sectional area, itr-^ 2 , for the purpose of computing drag 
and heat transfer to the body. The flatter the face, the more nearly constant 
are the variables, so the results would be most applicable to a flat -face body 
of area jtr-^ 2 . The application to a spherical body may be made in an approxi- 
mate way by introducing a shape factor, as described in the next section. The 
one -dimensional approach described above is equivalent to the approach used by 
Allen and James (ref. 3), who obtained consistently good results in general. 
However, Allen and James did not need to use the shape factor for the sphere, 
since they assumed an effective drag coefficient equal to unity, which is 
appropriate to a sphere, so that a shape factor was effectively "built in." 

In the present analysis, the drag coefficient will be calculated rather than 
assumed, so, as will be seen, the shape factor is needed here and its use 
yields consistent results. 

With this background, a procedure can now be formulated to calculate 
approximately the effects of the boundary shock wave on the meteor flight 
parameters. The altitude -density function, p, of reference 3 (see also 
table I) will be used, where 


p oo - p SL p 

is the ambient density at altitude h and 


(1) 
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( 2 ) 


P SL = 1 * 225 x 10 ~ 3 gm/cm 3 

is the ARDC standard atmosphere density at sea level (ref. 13 )* The density 
ratio across a normal shock -wave (needed in the calculations that follow) is 
approximately equal to the ratio of the ambient -stream density to the stagna- 
tion density behind the shock and, hence, is given for V > 7 . 5 X 10 5 cm/sec by 
the approximate relation (cf. ref. l 4 , p. 50) 


k 



-0.039 


v y ,X6 

10 2 cm/sec J 
__ - - 


( 3 ) 


To compute the heat radiation incident on the meteor face, one needs to 
know the volume of the outer shock layer. For a spherical nose at high hyper- 
sonic Mach number, the ratio of the shock standoff distance to the sphere 
radius is given approximately by Lighthill 1 s solution (ref. 15) • In the 
present case, the air -vapor interface would be the equivalent body, so that, 
assuming the Interface to be nearly spherical (see sketch (a)), one can 
approximate the outer standoff distance d Q by a form equivalent to 
Lighthill* s solution. Using the definition 

d^ 

— s D - 1 ( 4 ) 

r i 


where 


r 


i 


r b 


(5) 


one can obtain D implicitly as the appropriate root of 

A = 3(1 - k) 2 + 5 ( 4 k - 1)D 2 + 2(1 - k)(l - 6k) D 5 = 0 


(6) 


For the radiation rate per unit volume of the outer shock layer, which we 
denote as I, the relation given in reference lL, approximately /alid for 
V > I3.7XIO 5 cm/sec, is 

5.05 

U ^ -rt = 6 . 411X10 - 5 p 1 - 8 ( r- a V 1 (7) 

\1 erg cm 3 sec ly H ylO 2 cm/ secy 1 ' 


(Results obtained from use of this equation should be considered in the light 
of the qualifications made by the authors of reference l 4 .) In view of the 
discussion above, and with the assumption that one -half the radiation emitted 
by the outer shock layer impinges on the body, we represent the total radia- 
tive heating rate to a spherical body with a flattened face as 


1 p 

Q r = - otr b ^Id c 


( 8 ) 
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and define the radiative -heat -transfer coefficient 4 as 


/ \ Qr 

1 

T II X -rr3/ 2\ 

2 B J ^ ) 


Ido 
p V s 


( 9 ) 


Since d 0 /rp I s known from equations (4), (5), and ( 6 ), the value of 
(Cn /r^ )f>£ at any altitude for a flat -face tody is found from 


Ho 




ff 


p V 3 


vs- (D 


i) 


( 10 ) 


where V is to he obtained from the meteor data. 

The effect of the forces (other than body forces) that produce the 
changes in motion of the meteor may be represented parametrically by an effec- 
tive drag coefficient (cf. ref. 3),, defined by the relation 


C D ef f 2 Pco V2 (^b 2 ) = 


m 


dt 


= m - 


dV 

dt 


+ g sin Q 


( 11 ) 


where the meteor mass m is approximately 


4 3 

m = — rtr-fo p m 


(12) 


cm 


and p m is the mass density of the meteor. (See table II for the physical 
constants needed in this and other sections below. The value p m = 3.5 gm/c 
corresponds to the Ondrejov meteor Pribram (ref. 3) and is also nearly equal 
to the value 3.4 used by Opik (ref. 12) for typical stone meteors • ) From 
equations (ll) and (l2), Cp e ^ f /r b is known: 


"b e ff 8 p : 


m 


r b 


3 Pj 2 


dV 

dt 


+ g sin 6 


( 13 ) 


In the calculation it is convenient to use the quantity 


z = 


vc Hr 




c D e ff Cp 


( 14 ) 


'eff 




4 The reflectivity of the vaporizing surface material was disregarded 
because it is not known. 
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where Cr t ? for a "body that does not have a flat face. Is related to (^H r )ff • 
That relationship will he considered later. An expression for p e , the 
pressure outside the boundary shock wave, will also he needed. Note first 
that the pressure at the stagnation point 0 is the stagnation pressure along 
the line 2-0-e in sketch (h) and is given approximately hy 

p 0 = PcJ 2 (15) 


Then the ratio of the pressure p e to the stagnation pressure is given in 
terms of the Mach number outside the boundary shock wave. Mg, and the ratio of 
specific heats in the vapor, 7 (assumed constant), as (e.g., cf. ref. 16, 

p. 53 ) 


Pe 

Pq 




M c 


-7 


( 16 ) 


Based on the general ideas and formulations developed in this section, 
the equations needed to determine the flight parameters will be presented in 
the following sections, which include a summary of the calculation procedure. 


Shape Factor 

If the one -dimensional analysis with constant variables on the meteor 
face is used for flow near the meteor surface, the resulting equations are 
most appropriate to a flat -face body unless an assumed value of C B e ff 
appropriate to a sphere is used, as done by Allen and James in reference 3. 
When the value of ^R e ff n °t assumed, an appropriate shape factor must be 

used, if the equations are to be applied to a nearly spherical body, to insure 
realistic results . 

To derive a reasonably realistic shape factor, consider the following: 
Suppose a sphere and a flat -face body have the same cross-sectional area, ?rr 2 , 
and that some variable \| / on the surface has a constant value on the 

flat face normal to the flow and varies as 

V = i Q COS n cp (17) 

over the front half of the sphere, where cp is the angle between the axis, 
or stream direction, and a radial line of the sphere through a point on the 
surface. Let if be the integrated value of \|r over the front surface area. 
Then, for the flat face. 




ff 


^+0 


( 18 ) 


and, for the sphere. 
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¥ 


sphere 


n prt /2 

/ t dS = / cos n cp(2rtr 2 sin cp dcp) 

S o 


^TC/2 


= jtr 2 ^ 


2 cos 21 cp sin 9 dcp 




n 






(19) 


Suppose now one wants to compute the total radiative heat transfer to the 
surface for a heat flux that varies approximately as \|r Q cos n cp. The 
radiative -heat -transfer coefficients would then he related by 


2 

(^H r ) sphere n + 1 


(20a) 


where (Ojj r )^^, is given by equation ( 9 ) or (10). Now let the radiative -heat - 
transfer coefficient to he calculated for the meteoric body be represented by 


C Hr = ^ (C Hr ) ff (20b) 

where sjj is a shape factor for heat transfer. If, for example, the body is 
a sphere and if the radiative heating rate varies as cos 3 cp, for which n = 3 
in equation (20a) (cf. fig. 9 of ref. 17) * then Sjj = 2. Thus, 


s H = 2 for sphere with q r /q ro = cos 3 cp 
Sjj = 1 for flat -face body 


(21a) 


Similarly, suppose one wants to compute drag due to a force per unit area 
normal to the surface that varies approximately as \|r Q cos 2 cp (a reasonable 
assumption; true for Newtonian flow) . Then the force per unit area parallel 
to the axis would be ^ cos 3 cp. Using equation (l9)> one then finds the 

drag of a sphere to be one -half the value on the flat face (since n = 3) . 
Thus, assume 


CD eff 



efPff 


where the shape factor for drag, s^, is 

s-p = 2 for sphere 

s-q = 1 for flat -face body 


( 22 ) 


(21b) 
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Further, for convenience and in view of equations ( 21 a) and ( 21 b), consider a 
common shape factor for heat transfer and drag: 


S = Sr = S D 


(21c) 


An appropriate value for s may he determined as that value which causes 
C^eff *^ e un: *^y at altitude. For the Pribram meteor, 

s = 2 (21d) 


gives this result. 

The quantity Z 


in equation (l4) 


Z = 


VCH r 

CD eff 


is now 



Cr) eff 




(IV) 


Conservation Equations for Meteor Motion and Heating 

Now consider an equivalent flat -face body with uniform values of flow 
parameters over the meteor face, for which we want to write one -dimensional 
flow conservation equations. The one -dimensional equations of conservation 
of mass, momentum, and energy of the meteoric body can be developed in several 
ways. One derivation is given in appendix B for conditions and assumptions 
corresponding to the description in the above section. In appendix B the 
conservation principles are combined into equations (B8a) and (B 28 ) which may 
be designated, respectively, as the momentum equation and the energy equation. 
These may be written, neglecting P^/p s *- n comparison to unity, as 



where £ is given by equation (Bl6) as 


£ab = c sol( T f _ T a) + L f + c liq.( T b " T f) + L v ( 2 5) 


where T a is the cold interior temperature of the meteor. The cross- 
sectional area A in equation (23) Is itr-^ 2 , in line with the above discus- 
sion. The subscript b denotes the value of a quantity in the vapor flow at 
the boundary, the molten surface of the meteor. The notation is defined in 
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appendix A, but we point out here for convenience that u-^ is the vapor 
velocity relative to the boundary surface, the viscous compressive stress 

at the boundary, and q^ the associated heat conduction at the boundary 
(see ref. l) . (Some of the previous equations are discussed in appendix B; 
further results are discussed in appendix C.) 


To put these equations into a form that can be used in the present 
calculation, we use the definitions (cf. eqs. (9) and (ll)): 


(CH r ) ff 



( c D e ff) f f - 



1 

2 


P/A 


Jff 


(2 6 ) 


It is also convenient to make the approximation 


q c = -m 

and, in particular, 

^e-b = T b u b 

Equations (27) are exactly true for 

PC-n V.3 ^ 3 J 


k 


\ 1 \ 
- p + r K 


= 1 


(27a) 

(27b) 


(28) 


and are approximately true in general for Pr / 1 (see ref. l) , where p is 
the shear-viscosity coefficient, n the bulk viscosity coefficient, k the 
coefficient of thermal conductivity in the vapor, an d Cp the specific heat 
at constant pressure in the vapor. Note that if Pr = 1 and k = 0, the 
Prandtl number, pcp/k, is then 3/4. The vapor from a typical stone meteor 
may be approximated for some purposes as a perfect diatomic gas (see ref. 1), 
so that, in fact, Pr « 3/4 and Pr « 1 for this case. Equations (23) and (24-), 
substituted into (26) with use of (27b), take the forms: 


( Cl) efPff 



1 

7Mb 2 + 



(29) 


where 


(CH r ) ff 


g P > b u b 3 /lab l\ 
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( c H c ) ff 


(30) 


(31) 
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( 32 ) 
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Mb 2 » 


u-b 


7Pb/Pb 


( 33 ) 


and 7 is the ratio of specific heats in the vapor, Cp/c v . We may now obtain 
a second expression for Z (cf. eq. (l 4 )) by dividing VCn r from equation (30) 
by from (29) and ( 20 ). For convenience in performing an iterative 

calculation (to be described later), let us denote this expression by Z 1 and 
use equation (33) and the perfect gas equation of state 


Pb Pb RT b 

to eliminate u^, p^, and . The result is 


(3^) 


z* = 


VC Ht 


Mb 47 ^ (|~ + \ 7 M b ‘ 


CD 


eff 


1 + 7 Mb 2 + k 7Mb 2 C h 


( 35 ) 


which is to be set equal to 
calculation. 


Z in equation (lV) (p. 13 ) for the final 


Equations for Ablation Process and Boundary Shock Wave 

Abl ation process .- The temperature T^ at which the vaporization 
(boiling) takes place depends directly on the pressure to be overcome by the 
molecules escaping from the liquid (e.g., ref. l8, p. 88) and also on any 
significant final viscous stress (see appendix B of ref. l) ; the mass flux, or 
the rate of vapor ablation, depends directly on the number of liquid molecules 
near the surface that have sufficient energy to escape the liquid and on 
their velocities. The dependence of the temperature of vaporization T-^ on 
the total compressive stress, p^ - t^, for the phase transition at constant 
temperature is taken to be supplied approximately by the modified Clapeyron 1 s 
equation (appendix B in ref. l) in the form 

Ly / 1 1\ 

Pb - T b R v ref V 

— = e (36a) 

■^ref 

For the model of typical meteoric stone used by Opik (ref. 12 ), the mean 
normal boiling point (boiling temperature at a pressure of 760 mm Hg, or 
I.OI3XIO 6 dynes/cm 2 ) is 2960 ° K. This can be used as the reference point in 
equation (36a), so that, using L v and R from table II, we obtain the 
following expression for T^ as a function of p^ : 
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1° K 


36, 500 


26.13 Zn f ± dyne / cm 2 


' T b 


(36b) 


The relation used by Opik (ref. 12, p. 24) to represent the vaporization 
rate is 


(pu) b = | p b v x 


( 37 a) 


where (l/2)v x is the average molecular-velocity component normal to the 
surface in the +x direction. For meteoric stone (ref. 12, p. l6l), this 
reduces to 


^b 


1 gm cm“ 2 sec x 


= 3-08x10 


-4 


Pb 


1 dyne cm 


-2 


1 ° K 
r b / 


1/2 


(37b) 


which, with use of equation (34), becomes simply 

1/2. 




u b fR f Pb u b"b 


V/RTb 7 \ P i 


“ 3-08X10- 4 ' °' 333 


(38) 


where the constants from table IX have been used. 

Boundary sho ck wave . - At this point we note that, with equations (36b) 
and (38), Z 1 in equation (35) is determined in terms of p^ and Ch c j which 
are still unknown. Since is known and p e is known in terms of Mq 

(eq. (l6)), it is obvious that equations for the conditions across the bound- 
ary shock wave (i.e., conditions relating Me.? P e j etc., to etc.) 

and for Ch c are needed. If we define the density ratio across the boundary 
shock wave as 


£e _ 1 

^b " a 

we can write the following relationships (developed in ref. 1) : 



(39) 


(40) 
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(41) 


^- 1 + 3L^VU-«?) 

Pe _ ^ 

Pb a M e 2 


(42) 


Ch c 


2 




(43) 


Construction of Data Curve 

To use the equations given above for calculating the meteor flight 
parameters, one needs to know the values for velocity, V, and acceleration, 
dV/dt, as well as g sin 0, at each altitude, h. One of the most difficult 
problems of a study such as the present one is how best to use the meteor - 
tracking data. The data for velocity versus altitude are given in table III 
(taken from ref. 4, table 7). The data curve (fitted to the data points) 
must give physically reasonable results. Least-squares fits (see ref. 3) have 
given realistic results in the past except in the case of Ondrejov meteor 
Pribram, for which there was some difficulty. The shape of the curve fitted 
to the data appears to be especially important in the upper part of the tra- 
jectory, where the accelerations are small. If one data point is signifi- 
cantly in error or if some supposedly small factor is neglected, the results 
can be affected significantly. 

The view is taken here that the data curve must best fit the tracking 
data and at the same time yield a physically possible result (in particular, 
a nonincreasing radius) . In constructing the data curves of h versus V, 
certain appropriate assumptions can be made that allow one to calculate 
approximately the slopes near the top of the trajectory for different values 
of an assumed initial radius, r^*. Then, using these slopes as guides, one 
can graphically fit a curve to the data. This is, in effect, a combination of 
fitting a curve graphically to the data and using "approximate isoclines” to 
determine the shape of a portion of the curve. Of the resulting curves for 
different r^*, that value is used for which the results are most consistent 
along with the best apparent fit to the data. 

The procedure for constructing the data curve is as follows: The data, 

given in table III, are plotted in figures l(a) through (f). Note from 
table III that V does not change significantly, especially between, say, 
h = 90 and 6 5 km. For this altitude range, assume V = V* = 20.87X10 5 cm/sec. 
Also assume ^D e ff <3-°es not change much from unity and does not change 

much from its initial value, r^*, in this altitude range. These assumptions 
will be justified by the calculated results. From reference 3, sin 0 = 0.6853- 
From equation (13) we then have the approximate relation for the top part of 
the trajectory (90 to 65 km altitude) : 
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I km sec -1 


I km sec' 




* = 35 cm ( f ) r b * = 50 0 

Figure 1. - Graphical fit to tracking data of 
Ondrejov meteor Pribram. 



Using also 
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dV 

dt 


-V sin 9 

an 


(45) 


we then obtain for the approximate slope of the h, V curves: 


or 


dV = 1 _ A ^SL pV* 2 
dh V* \8 P ra r b * sin 9 
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h 


-lkm- 


< 90 


(46a) 


dV/dh _ 399 p 

1 sec (r b */i cm) 


4.70X10 4 , 


(65 < ^ < 9o) 


(46b) 


With p (from table I) and various values of r ^* , approximate values of 
dV/dh were calculated from equation (46b) and are listed in table IV. From 
these values, slopes were drawn on parts (a) through (f) of figure 1. The 
best curve was drawn through the points guided by these slopes. Since two 
data points (data from two different cameras) are given for h = 76 km that 
are very close, those points were assumed to be accurate. On the other hand, 
the top data point (h = 88 km) was virtually ignored. It was impossible to 
obtain realistic results using the top data point. 5 The data to be used as 
input (V and dV/dt) were then obtained from figure 1 for each r^* by meas- 
uring the slopes and using equation (45) - The values are listed in table V. 


Calculation of Flight Parameters 

For convenience, sketch (c) shows the overall scheme for calculating the 
flight parameters. 

Some parameters of interest that can be calculated, in addition to those 
already defined, are: the pressure drag coefficient, 

C Dp = (8/2^2 (4 7 

5 Wote that the data curve of Allen and James (ref. 3 for the same data) 
also appears to disregard the top point. Note also that the two top data 
points of table III were obtained with a camera different from that used to 
obtain the remaining points. 
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Meteor -tracking data: h versus V and dV/dt (curve graphically 

fitted to the data) ; and sin Q 

Altitude -density function: p versus h 

Approximations for: 

Shock density ratio, k 
Shock standoff function, D 

Radiation rate per unit volume of hot air layer, I 
Stagnation pressure, p Q 

Isentropic relation between stagnation pressure, p Q , and pressure 
outside the boundary shock wave, p e 

One -dimensional conservation equations for meteor motion and 
heating, using the shape factor, s, for sphere 

Vaporization equations and stone -meteor constants from Opik*s data 


Equations for conditions across the boundary shock wave 



Output: Altitude versus p^ T b , u b , P e Av P e /Pb> T e/ T b; M e> 

C-H r ^ ^Deff'* > ODp; and Cp^_ 


Sketch (c).- Summary of calculation. 

the drag coefficient due to the viscous -compressive stress. 


“Ti 


b v r 1 v ^H c ) ff 

Cd vc s(l/2) Poo V 2 s(l/2)p V® su b Lri/2) Poo V3j SUb 




the contribution to the total ^D e ff due to the retrothrust of the vapor, 

2 


Pb u b 


Cri t s(l/2)p V 2 Cr) eff Cl) p " Cd vc 


(49) 


and the heat -transfer coefficients: 


C H c s 


Ch - CH r + C Hc 


(50) 

(5D 


In equations ( 47 ) through (50), the shape factor has been used for consistency 
with equations (20b) and (22). 
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The detailed calculation procedure, including several simple auxiliary 
expressions that need no comment, is outlined in appendix D. The calculations 
were performed on an IBM 709^- electronic data processing machine. 


RESULTS AND DISCUSSION 


As outlined under "Construction of Data Curve," the Initial meteor radius 
is that value of the assumed initial radius, which gives the most consis- 

tent results along with the best data fit. The data curves for various r^* 
are given in figure 1, and the resulting curves for variations of r^ with h 
for the various r^* (calculated as outlined above and in appendix D) are 
shown in figure 2. For r^* less than 20 cm, the radius appeared to increase 
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Figure 2.- Variation of meteor radius with altitude (s = 2). 
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before decreasing. For r^* greater than 25 cm, the data fit does not appear 
to be as good as for 20 cm. It appears that r^* = 20 cm gives the best 
results, so the remaining results given here correspond to = 20 cm. (The 

data, listed in table V(c), are taken from figure 1(c).) Although the size of 
the main body of the meteor, before its fragmentation at lower altitudes, is 
very uncertain (see ref. 4), it is believed that the radius computed here is 
too low. This reflects, certainly, the approximate nature of these calcula- 
tions and, probably, inaccuracy and uncertainties in the original meteor data. 
The results here would be improved if the shape factor, s, decreased from two 
toward unity with decreasing altitude, so that would increase; 

would then decrease more slowly and have a larger final value. 


h 

I km 


\ 


The results indicate that if the idealizations in this calculation are 
approximately valid, then the rapid variations of flow variables and high heat 
conduction characterizing a strong boundary shock wave were present in the 
ablating vapor flow of the Pribram meteor. This conclusion is indicated by: 

(l) the heat -conduction parameter, (see fig. 3), is different from zero, 

(2) the ratios of conditions across the 
boundary shock wave (e.g., p e /p^; see 
fig. k) are greater than unity, and 

(3) the vapor -flow Reynolds number, 
Pb^-b^i/Pb (where <*i and were 
estimated roughly; see appendix E), is 
much greater than unity. (For alti- 
tudes less than about 65 km, 
is of order 10 2 or greater; see 
appendix E. ) The boundary shock wave 
is considered to be strong because of 
the large values of Cfr and Pe/p^. 

The calculated temperature change 
across the boundary shock wave on the 
Pribram meteor is small because the 
vapor ablates at low subsonic Mach 
number (M^ = 1/3) and because the den- 
sity ratio across the boundary shock 
wave is high (see eq. (tl)). 



v 

w 



2 3 

Temperature, velocity, heat -conduction parameter 

(a) Temperature, velocity, heat -conduct ion 
parameter . 



(b) Pressure, density. 
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Figure 3.- Vapor-flow variables at the meteor- 
vapor boundary (r^* = 20 cm, s = 2) . 


Figure 4.- Conditions across the boundary shock 
wave on the meteor (r^* = 20 cm, s = 2) . 


The validity of the equations used for the boundary shock wave depends on 
Re^ = p^u-^d^/fl^ being large (» l) . The results that depend on the boundary- 
shock -wave theory are therefore not reliable much above 6 5 km altitude in 
these calculations for the Pribram meteor. In the calculations below 6 5 km 
altitude, for which Rep is large, the vapor flow necessarily divides into 
thin viscous regions (the boundary shock wave and the interfacial layer) and 
an essentially inviscid region (see sketch (b )). 6 All viscous effects, 
including convective heat transfer, are absent from the inviscid region. Con- 
vective heat transfer from the hot compressed air is therefore confined to a 
thin layer at the interface; in the usual terminology of ablation theory, the 
convective heating is completely "blocked” from the body by the efflux of 
vapor. The only heat transfer from the outer layer of hot air behind the 
shock wave to reach the body is assumed to be that due to radiation. However, 
there is significant heat conduction in the vapor at the surface (within the 
thin boundary shock wave) that will be discussed below. Just as the boundary- 
shock-wave theory is not applicable at the beginning of the entry into the 
atmosphere, some point may be reached later where it is again not applicable. 

A point may be reached where the radius rp has decreased to such an extent 
that (l) Rep is decreased signif icantly because it contains a length propor- 
tional to rp and ( 2 ) the radiation is much lower for smaller rp so that 
the mass flux ppUp contained in Rep is also decreased significantly. 

Hence, when the radius becomes small, the vapor flow near the body may again 
approach a boundary-layer -type flow. This point in the trajectory of the 
Pribram meteor is below the lowest data point (see appendix E for values of 
Rep), if it occurred at all. An additional effect of this phenomenon, when 
the radius becomes small enough and/or the velocity slow enough so that the 
ablation rate decreases considerably, is that the radius will decrease much 
slower; that is, the curve of radius versus altitude will at some point again 
become more nearly vertical in figure 2. This may be the effect beginning to 
show up in the curve for rp* = 25 cm in figure 2 . 

The primary meteor flight parameters of interest are the radius (fig. 2), 
the drag coefficients, and the heat -transfer coefficients (see figs. 5 ,? 6 , 
and 7 ). The first important result regarding drag to be noticed in figure 5 
is that c D e pp remains, for all practical purposes, unity. (This depends, of 

course, on the shape factor being constant at 2.) Note, then, that one could 
ignore the boundary shock wave, assume = 1 . 0 , and determine rp from 

equation ( 13 ) (following the general procedure of Allen and James, ref. 3 ), 
using the most appropriate data curve for V and dV/dt versus altitude. The 
radius and Cp ^ would then be known accurately, but not necessarily the 
heat transfer (see below) or the thermodynamic properties of the vapor. The 
second result regarding drag is that the viscous -compressive stress (at the 
vapor boundary) produces most of the effective drag, the thermodynamic 

6 Note that, even at 65 km altitude (taken to be the upper limit of 
validity of the theory), where the meteor radius was about 20 cm (fig. 2 ) and 
the thickness of the vapor layer was about 0.0l6 X 20 cm = 0. 32 cm (see 
appendix E) , the ratio of the boundary shock thickness to the entire vapor- 
layer thickness (ratio of order l/Rep) was about l/77> so the thickness of the 
boundary shock wave was about O .32 cm X 1/77 = 0.00^15 cm. 
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Drag coefficients 

Figure 5 • - Meteor drag coefficients (r b * = 20 cm, 
s- 2 ). 
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Figure 6 .- Meteor heat -transfer coefficients 
(rb* = 20 cm, s = 2 ) . 



Comparison parameters 


Figure 7 .- Ratios showing boundary-shock effects 
(r b * = 20 cm, s - 2 ) . 

conduction in the boundary shock wave 
fig. 7 ) . To aid in understanding the 
at the surface, it is useful to consic 
have (see ref. l) : 


pressure drag being generally much 
smaller and the retrothrust (repre- 
sented by a M drag coefficient,” Cp^) 

practically negligible (see figs. 5 
and 7 ). This result, that Cp vc is 

the most significant portion of Cp ^ 
can be misleading, however. If the 
viscous-compressive stress were not 
included in the analysis, the pressure 
drag calculated would be much larger, 
since p^ would take the place of 
lb + (~ T b) in momentum equation 
(eq. (23)). 

The plot of the heat -transfer 
coefficients due to radiation and heat 
conduction in the vapor at the boundary 
(fig.. 6 ) indicates that the heat 
plays a substantial role (see also 
roles of radiation and heat conduction 
.er equation (24) . Assuming Pr = 1, we 


so that, from equation (24), 


% = T b u b 


^r 





C52) 


(53) 


2h 



From equations (53)> (B15)j and (Bl6) (or eqs. ( 74 ) in ref. 1), one sees that 
the radiative flux to the body, -q r , contributes to: (l) heating the body 

material to vaporization temperature (by conduction within the body), ( 2 ) 
vaporizing the body material, and ( 3 ) increasing the kinetic energy of the 
vaporized material [(1/2)^^] . Included in the energy of vaporization is the 
"flow energy," p^/p^ = RT^, added to the vaporized material. 'The latter term 

^Pb u b(lb/Pb) = u bPfc)J is equivalent to the work required to blow off the vapor 
against the pressure p^ with velocity u^ • This energy, added by work, 
stays with the vapor as flow energy (giving the vapor an "enthalpy," e + p/p), 
since the force due to the thermodynamic pressure p is not dissipative. Now 
we see from equation ( 52 ), by analogy, that the contribution of the heat con- 
duction back through the gas at the boundary, -q^, is in the form of the work 
required to blow off the vapor against the viscous -compressive stress, -t-^, 
with velocity u^. Since the force due to -t^ is dissipative, however, no 
energy is added to the gas by this work, and it does not appear in 
equation ( 53 ) • 

For the heating, then, although the viscous effect (the heat conduction 
in the boundary shock wave) is large, one sees (as in the case of drag) that 
omitting the boundary shock wave from the analysis would not influence the 
overall results as much as might at first be supposed. The heat conduction 
simply balances the work of blowing off the vapor against the viscous - 
compressive stress (assuming Pr = l) . 

Although detailed flow calculations would be affected significantly by 
the boundary shock wave, an important result here is that, for overall gross 
effects, ignoring the boundary shock wave and assuming CD e ff to be unity and 
the heat load producing the vaporization to be just the radiation from the hot 
gas, one should obtain good approximate results. The only difficulty remain- 
ing in ignoring the boundary shock wave is in determining the ablation temper- 
ature, which also determines . However, since the largest part of 

is comprised of the heat of vaporization, is not too strongly dependent 

on the vaporization temperature. The approximate equations used by Allen and 
James (ref. 3) therefore give good approximate results for the overall motion 
and heating of a meteor (especially since their assumed values of £ and T-^ 
are approximately correct for the most critical part of the trajectory), even 
though the boundary shock wave has not previously been considered. The most 
important conclusion here is that detailed flow calculations, both for 
determining the flow characteristics and for computing more precise results on 
the motion and heating of the meteoric body, should consider the possible 
effects of the boundary shock wave. In the case of the Pribram meteor, the 
boundary shock effects on the flow are indicated by the idealized theory to be 
substantial. 


CONCLUDING REMARKS 


The results of including a boundary shock wave in a simplified 
calculation of the flight parameters of a stone-meteoric fireball have been 
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studied using meteor-tracking data for the analysis. The results found or 
deduced are enumerated explicitly, and then elaborated upon, in summary of the 
above discussion. They include: 

(a) the calculated values of ratios (which represent the "jumps" or 
rapid variations ) of vapor-flow parameters across the boundary 
shock wave, and values of other parameters in the vapor that 
relate to the boundary shock wave; 

(b) qualitative differences in values of the vapor-flow parameters 
from values that would be obtained by ignoring the possible 
presence of the boundary shock wave (these differences can be 
inferred from the results (a), but were also checked in a simple 
calculation, not shown); 

(c) the significance of (a) and (b) to the vapor flow field; 

(d) the calculated values of the flight parameters, including the 
meteor radius and the various heat-transfer coefficients and 
drag coefficients; 

(e) qualitative differences in the calculated values of the overall 
flight parameters from values that would be obtained by ignoring 
the possible presence of the boundary shock wave; and 

(f ) significance of the above results (of including the boundary 
shock wave in the present calculation of meteor flight parameters) 
to more detailed analyses of the flight parameters. 

The calculations indicate that, if the idealizations made are valid, 
conditions were present in the ablating -vapor flow of the Ondrejov meteor 
Pribram for a strong boundary shock to occur according to the boundary -shock- 
wave theory presented previously. 

V ' 

During the portion of the Pribram meteor trajectory for which tracking 
data were obtained and for which the boundary- shock- wave theory is presumed to 
apply (between about 65 and 45 km altitude) the results are: (l) neither the 

pressure ratio nor the density ratio across the boundary shock wave was ever 
less than about 4.5, ( 2 ) the temperature change was about 2 percent across the 
boundary shock, ( 3 ) the viscous -compressive drag was always at least 75 per- 
cent of the total effective drag, (4) the total effective drag coefficient 
remained nearly constant at the value unity, and ( 5 ) heat transfer by conduc- 
tion at the surface in the vapor was always at least 20 percent as large as 
the radiative heat transfer. 

Since the density a short distance from the molten surface, p e , is 
significantly different from p^, which in turn is fairly close to the density 
that would be calculated by ignoring the possible presence of a boundary shock 
(or equivalently by omitting the translational-nonequilibrium condition on the 
vaporization rate), the calculated effects on the vapor flow are substantial • 
Therefore a more detailed analysis of the flight parameters that would use a 
detailed flow field analysis (for more precise calculation or for other 
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II 


reasons) should include consideration of a boundary shock vave ; under the 
conditions where its presence in significant strength is expected. 

The results of the analysis including the boundary shock wave indicate an 
initial radius between 20 and 25 cm ; probably closer to 20 cm, for the meteor. 
The actual initial radius is believed to have been somewhat larger than this. 
The difference is probably due to the approximations made in the analysis and 
to uncertainties in the original meteor data. 

If the boundary shock wave were not considered in calculating the overall 
gross flight parameters (so that virtually all of the heating would be assumed 
to be by radiation and most of the drag due to pressure), it is still possible 
to obtain realistic approximate results for the overall motion and heating of 
the meteoric body by judicious choice of: (l) assumed values of the effective 

drag coefficient (CD e ff)> the temperature of vaporization (T^), and the total 
energy required to heat and vaporize a unit mass of the meteoric material 
( £ a *b 7 (as used effectively by Allen and James) and (2) a suitable curve 
representing the meteor-tracking data. Thus, although the occurrence of a 
boundary shock wave may be important to the flow details, the overall flight 
characteristics of the meteor (in the class expected to have a boundary shock) 
are not too significantly affected by the presence of a boundary shock wave 
if proper estimates are made, as done by Allen and James. 
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APPENDIX A 
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PRINCIPAL NOTATION 
cross-sectional area of the body 

effective drag coefficient; equations (ll), (22) , and (26) 

coefficient of drag due to pressure, equation (k-j) 

contribution to total due netrothrust of vapor, 

equation ( 49 ) 

coefficient of drag due to viscous compressive stress, 
equation (48) 

coefficient of heat transfer by conduction, equations (32) 
and (50) 

conduction coefficient defined by equation (31) 

coefficient of heat transfer by radiation; equations (9), (20b), 
and (26) 

specific heat 

specific heats in gas (vapor) at constant pressure and constant 
volume 

shock standoff function defined by equation (4) 

distance from body to interface between vapor and air 

distance from interface to shock wave, that is, thickness of 
outer shock layer 

differential quantities first used In equation (BIO) 
differential quantities defined prior to equation (Bl) 
specific internal energy 

total force per unit area acting on an elemental surface 

body force per unit mass acting in the direction of motion; 
f = g sin 6 when 9 = gh 

acceleration of gravity 


h 

I 


Ki,K 2 
KiS K 2 « 


k 

k 

L f 

Lv 

M 

m 

m 

1TL 

P 

Pr 

Pr 

Qr 

q. 

% 

R 

Re b 

r b* 


altitude above the earth’s surface 

rate of energy radiation per unit volume in the outer shock layer 

constants defined by equations (C5) or (C6) 

constants of integration first used in equations (C9) and. 
evaluated in (Cll) or (C12) 

thermal conductivity 

shock density ratio, P^/p^ 

latent heat of fusion 

latent heat for vaporization at low rate 
Mach number 
mass of meteor 

portion of body mass between x a and x = 0 at a given time 

molecular weight of vapor 

pressure 

Prandtl number, pcp/k 

Prandtl number based on p, pcp/k 

total radiative heating rate to a flat -face body 

heat flux (positive to the left in sketches (a) and (b), positive 
to the right in sketches (d) and (e ) ) 

conduction heat flux 

radiative heat flux 

gas constant (universal gas constant divided by molecular weight) 

vapor -ablation Reynolds number, Pb u bd-i/£t> 

body radius, sketch (a) 

assumed initial radius of meteor 

distance from center of nearly spherical body to vapor -air 
interface , sketch ( a) 


s 


shape factor, equations (21) 


T temperature 

t time 

u velocity of vapor relative to molten meteor surface (to the right in 

sketch (d)) 

V magnitude of meteor velocity relative to the earth 

V* approximate value of meteor velocity at atmosphere entry 

Vfc velocity of vapor at the molten surface relative to the body mass 

whose velocity is V; u^ - u g 

x coordinate in the direction of the body motion whose origin is fixed 

at the receding molten surface, sketch (d) 

x 1 coordinate in the direction of the body motion whose origin is fixed 

in space 


Z,Z T , equations (l 4 ) and (35) 

UD eff 


a 

PijPa 

7 

^ab 

0 

7\ 

P 

P 

Pm 

P SL 

P 


— , equation (39) 

p e 

constants defined by equations (C 4 ) 

c-p 

ratio of specific heats, — 

c v 

e b “ e a + equations (B 15 ) and (Bl6) 

angle between velocity vector and equipotential surface (e.g., 
earth* s surface when cp = gh) 

function of k and D, equation (6) 

shear-viscosity coefficient 

defined in equation (28) 

mass density per unit volume 

mass density of meteor 

atmospheric density at sea level, equation (2) 
altitude -density function, equation (l) 
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a Z - Z 1 

t viscous compressive stress 

cp body-force potential; for gravitational force cp = gh, equation (Btb) 

f variable quantity on meteor surface , equation (17) 

\|r Q constant value of \|r on a flat -face body 

Subscripts 

a in general, an arbitrary point in x < 0; in particular, value where 

dT 

k — & 0, that is, where T a is the cold interior temperature 
b value at the boundary in the vapor (x = 0 + ) 

e value outside the boundary shock wave in the vapor 

f value at fusion, that is, where T = Tf, the fusion temperature 

ff value on flat -face body 

liq value in the liquid (molten meteoric material) 

o value at the stagnation point, on the air-vapor interface 

ref arbitrary reference value 

s value at the molten surface, inside the surface (x = 0~) 

sol value in the solid meteoric material 

2 value behind the shock wave, sketch (b) 

00 condition in the ambient atmosphere at altitude h 
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APPENDIX B 


DERIVATION OF THE ONE -DIMENSIONAL CONSERVATION EQUATIONS 


A derivation of the equations of conservation of mass, momentum, and 
energy is presented for an accelerating body with the following conditions: A 

rapid efflux of vapor from the forward face of the body is produced by vapor- 
ization of the body material due to absorption at the body surface of intense 
radiation. All motion, including that of the body and the vapor, is assumed 
to be one -dimensional. The body has a cross-sectional area A and an 
instantaneous mass m. The gas in front of the body is the outflowing vapor. 
The effects of any gas to the rear are neglected. The body is moving in a 
force field of potential cp. (in the case of gravitational force, cp = gh. ) 

For the following development, refer to sketch (d), where all motion is 
to the right. Let 

V = V(t) be the instantaneous velocity of the body; 

x be a coordinate in the direction of the body motion whose 

origin is always fixed at the receding surface; 

x f be a coordinate in the direction of the body motion whose 

origin is fixed in space (e.g., with respect to the earth* s 
surface) ; 

u s = u s (t) be the rate at which the surface is receding into the body mass, 
that is, the velocity in the -x direction of the surface 
relative to the internal mass of the body. (Thus, the 
absolute velocity of the surface is V - u s . ) ; 

u-fr = Ub(t) be the velocity of the gas (vapor) at the surface relative to 

the surface. (The absolute instantaneous velocity of the gas 
at the surface is - u s + V.); and 

u = u(x,t) be the instantaneous velocity of the gas to the right relative 
to the surface at a distance x from the surface. (The 
absolute velocity of the gas at x,t is u - u s + V.). 

Note that, relative to the x coordinate, the surface is stationary, the mass 
inside the body is moving with velocity u s , the gas is moving with velocity 
u, and the surface value is u-^ . (Relative to the x* coordinate, the body 
has velocity V = dx f /dt.) (Subscript b refers to values in the gas at the 
boundary of the gas flow (x = 0 + ) ; subscript s refers to values inside the 
unvaporized material at the surface (x = 0”) ; subscript a refers to an 
arbitrary point a inside the surface (x < 0).) 

Consider the motion of the system of constant mass m during the 
infinitesimal time interval dt between t Q and t Q + dt, where m is the 
portion of the body mass m between x a and x = 0 at time t Q . The mass 


32 



Time t = t 0 


h 


-d€ 


... 


" 1 


Stress fb 

Gas velocity ub -u s + V 




Time 1 = t n + dt 


f« + df a -c 


Solid and liquid mass 
m with velocity V 


s' a b 


s 

p* ' 


Surface x r 0 moving 
with velocity V - u s 


+ dp b 


Gas velocity u b - u s + V + 
d(u b - u s + V) 


So 1 d and '-quid moss m+dm 
w th veiocify V + dV 


Surface * - 0 moving with 
velocity V-u s + d(V -u s ) 


Stress fb + df b + dij 


Velociiy u b -u s + V 

+ d (ub~u s + V) 

+ f^ d7? 


Sketch ( d) . - Ablation of surface during time interval dt. 


vaporized during dt is -dm. (Although -dm is lost from the "body" during 
dt, it remains part of our "closed, " or constant mass, system. ) Let be 
the dimension of the mass -dm before it is vaporized and dr] its dimension 
after it is vaporized; thus. 


where, by definition. 


-dm = p s A d£ = p^A d r\ 


u 


di 

dt > 


u b = 


dr) 

dt 


(Bl) 


(B2) 
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From equations (Bl) and (B2) we have a convenient statement of conservation of 
mass : 


“ dt = Ap s u s = A % u b (B3) 

In the following an appropriate "momentum equation" is derived from the 
principle of conservation of momentum. First, define f to he the sum of 
"surface forces" per unit area (such as pressure and viscous stress in the gas 
or liquid, and tensile or compressive stress in the solid) . Positive values of 
f indicate tension, so that, in the gas or liquid, 

f = “P + t (B4a) 

where p is the thermodynamic pressure and t the viscous -compressive stress 
(t = ji(du/dx); see, e.g., ref. 16, p. 33 I, or ref. l) . Also, define f to he 
the hody force per unit mass acting on the mass m in the direction of the 
hody motion, that is, 

1 “ ax' 



(in the case where the potential is cp = gh, 

f = g sin Q 


where 


sin G 


dh 

dx* 


dh 
V dt 


i.e., where G is the angle between the trajectory of the hody motion and an 
equipotential surface, such as the earth's surface.) The average external 
force to the right (sketch (d)) acting on the mass m during dt is then 


mf + ^ A 


(f b - f a ) + (f b + df b + jg- - (f a + df a ) 


= mf + A(f h - f a ) + higher order terms 


so that the impulse during 


The momentum 
(m + dm) (V + 


change (to the 


dV) + (-dm) 


U v 


dt is 

mf dt + A(f-j 3 - f a )dt 
right) of the mass m during dt is 

- u s + V + d(u b - u s + V) + | 1^- dT) 

= m dV - (u-^ - u s )dm + higher order 


- mV 
terms 


(B5) 


(B 6 ) 
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Neglecting the higher order different ials, we may then write, from the 
principle of conservation of momentum. 


A(f- b - f a )dt + mf dt = m dV - (u b - u s )dm 


(B7a) 


Convenient forms for this equation are 


m /dV 
A \dt 



f a - f h + ( u b " u s) 


1 dm\ 
A dtj 


(B7b) 


and, with use of equation (B 3 ), 



m dV 
A dt 



Ph u h 2 " f h 


(B7c) 


If we take point a as the hack surface of the "body, where f a = 0 and m = m, 
we have, using equations ( B 3 ) and (B4) in(B7b), 


m / dV 
A \dt 




(B 8 a) 


in which, of course, P^/p s « 1 and may be neglected. Another convenient 
form for equation (B 8 a) is 


-m 



Bp + B-yc + 


(B8b) 


where the pressure drag, the viscous -compressive drag, and the retro thrust due 
to ejected vapor are, respectively, 

D p = Pfc)A 


■^vc 


Fiji = 


-T b A 


dm / \ dm 

" (u b ' u s) = “ dt 


(B9) 


“ Ph% 2A ( 1 - 7T 


f>b u b 2A 


where Vp = Up - u g is the gas velocity relative to the body mass whose 
velocity is V. 


To derive an equation representing conservation of energy, we begin with 
the first law of thermodynamics applied to the closed system of mass m 
during dt (see above definitions and sketch (d) ) : 
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dE = dH + dW (BIO) 

where 

dE = the change in energy of the system of mass m during dt 

dH = the net heat added to the system of mass m during dt 

dW = the work done on the system of mass m hy external forces during dt 

The energy E possessed by the mass m is composed of internal energy, 
kinetic energy, and potential energy. Let e denote internal energy per unit 
mass at a certain point in the mass m at time t = t Q and e + de be the 
corresponding value at time t Q + dt. Also, denote the latent heat of fusion 
and latent heat of vaporization at low rates (see appendix B of ref. 1) by 

Lf and L v , respectively. Refer to sketch (e) . During dt the surface, x = 0, 



Sketch (e).- Change of specif ic -internal -energy distribution in body during time interval dt. 

where the vaporization is taking place, recedes from s to s T , a distance d£ 
relative to the internal body mass. For simplicity, changes in e s and e-^ 

with time are neglected, that is, e s r ^ e s and e^t e^ . Then, during dt, 

the internal energy curve is simply shifted a distance d£ to the left in 
sketch (e) . Neglecting higher order differentials, one finds the shaded area 
to be 

( e s - e a) d i 

Since the mass included between point a and point s f is m + dm, the 
change in internal energy of the mass m + dm during dt is 

P s A ( e s - e a ) d| = -dm(e s - e a ) (Bll) 

Note that this takes into account a phase change (fusion) at x = xf where 
the fusion temperature is Tf ; thus, if the specific heats in the solid and 

liquid are constants c so q and C]_j^, one may write 
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e s e a “ c sol^f + Lf + c liq^-'-s “ ^-f) 

At t = t Q , the mass (-dm) has an internal energy of 

-dm (e s “ ^ ^ + higher order terms 


(B12) 


and at t = t Q + dt the mass ( -dm) has an internal energy of 


-dm ( ~ — dr\ \ + higher order terms 


so that, with neglect of higher order differentials, the change in internal 
energy of the mass (-dm) is 

-dm(e b - e s ) (BI 3 ) 

Thus, from equations (Bll) and (BI 3 ), the total change in internal energy of 
the mass m during dt is 


-dmUab - RT b ) 


(B14) 


where 


?ab - ET b ~ % - e a ( B1 5) 

(see appendix B of ref. l) . For constant specific heat in the solid, liquid, 
and const ant -temperature phase changes (e.g., see ref. 19 )^ if point a is 
to the left of the fusion point, then 

Cab = c S ol(*f - T a ) + % + c liq (T b - T f ) + L v , (x a < x f ) (Bl6) 

In sketch (d), the kinetic energy of m at t Q is (l/2)mV 2 . At time 
t Q + dt the kinetic energy of the mass m = (m + dm) + (-dm) is 


| (m + dm) (V + 


dV)‘ 


!<■ 


•dm) 


Therefore, the total change in kinetic energy of the mass m during dt 
(omitting higher order differentials) is 


mV dV + 


(-cam) 





(B17) 


From expressions (BlU), 
change : 


(BI 7 ), and the expression for the potential -energy 
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m dcp = -mf cbc' = -mfV dt 


we then have 


dE = mV dV + ( -dm) 


£ab - HTh + | (u b - u s ) 2 + V(u b - u s ) 


- mVf dt 


The net heat added to the system of mass m during dt is 


clH = (<i a - ct b )A dt 

where q denotes heat flux to the right in sketch (d) and where 

<l a = -(k dT/dx) a 

% = <i r + q. Cb = q. r - (k dT/dx) b 


(B18) 


(B19) 


(B20) 


and -q is the radiative heat from the right that is absorbed at x = 0 and 
q Cb the heat flux by conduction in the vapor at the boundary (x = 0 + ) . 


The work done on the system of mass m during dt by the gas in front 
of it is the product of the average force during dt on the forward bounding 
surface times the distance through which that surface is moved. From 
sketch (d) one finds that the forward bounding surface of the mass m moves, 
during dt, from b to c, an actual distance of 


(u b - u s + V)dt + higher order terms 


Thus, the work done on m during dt by the gas in front is 


f^A^u^ - u s + V)dt + higher order terms (B2l) 

Similarly, the work done on m during dt by the mass behind m is 


. -f a AV dt + higher order terms (B22) 

and the total work done on mass m during dt, with neglect of the higher 
order terms, is then 


£W = (f b - f a )AVdt+ f b (u b - u s )A dt (B23) 

The first law of thermodynamics, equation (BIO), written in the form 


dH = 
dt 

and with use of the expressions (Bl8), 
expressing conservation of energy: 


dE _ dW 
dt dt 

(B19), and (B23), 


(B24) 

gives an equation 
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' % = I 


( SL - A 

1 

dm 

W V 

" A 

dt 


^ab “ ^-b p ^ u b u s) + ^ u b “ u s^ 


+ (f - f„)V - f b (u„ - u s ) 


(B25) 


Substitution of the momentum equation (B7b) into (B25) yields 

" f b (u b " u s) 


1 dm 

q a % ' a dt 


^ab " RT b + 2 (u b " U a) 2 


(B26) 


Another useful form obtainable from equation (B25) with substitution of (B 3 ), 
(B7c), and (B15) is 


'Vb 6 ^ + 2 


^b u b 3 + 


% 


u b f b 


p s u s e a 


+ 2 p sV 


+ 


■p m dV m f 

a A dt A 


(B27) 


If we now take point a far enough back from the body surface that q a is 
negligible and T a is the cold interior temperature, using equations (B 3 ) 
and (b 4), we may write (B2 6 ) in the following two convenient forms: 


-(q r + i Cb ) 


= % u b 


5 ab 


RTb + 2 u b = 


h 

Pc 


2 -| 


+ (Pb " T b )u b (l - ^ 


(B28) 


where, as in equation (B 8 a), Pb/Pg « 1 and may be neglected; or 


dH 

dt 



RT b 



+ (D p + D vc )v- b 


(B29) 


where dH/dt = -q^A = -(q r + q c ^)A; v b = u b " u s ; bh e S as velocity relative to 
the body mass whose velocity is V; D p and D vc are given by equation (B9), 
dm/dt is given by equation (B 3 ); and £ is given by equation (Bl 6 ), with 
T a as the cold interior temperature. 

Following is an interpretation of equation (B29) : The effects of the 

heat transfer dH during the time interval dt are: (a) to add internal 

energy to vaporize the mass ( -dm) and raise the temperature of an equivalent 
mass to the level that ( -dm) had before being vaporized (this accounts for the 
term -dm^^ - RT^); (b) to add kinetic energy to (-dm) relative to its 
original motion [(l/2)(-dm) v^ 2 ] ; and (c) to blow off vapor with force (Dp + D vc ) 

acting through distance v^dt relative to the body motion [(Dp + B vc )vj d dt]. 
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In this derivation of the equations of motion and heating, note 
specifically that the following are accounted for: 

(1) heat transfer into the body, as well as heat transfer used to 
vaporize the surface material; 

(2) work done to blow off vapor, as well as the kinetic, potential, and 
internal energy possessed by the outflowing vapor; 

(3) recession of the ablating surface; 

(4) drag due to the viscous compressive stress, as well as pressure drag, 
retrothrust of ejected vapor, and the gravitational force; 

(5) heat conduction associated with the viscous stress, as well as 
radiative heating. 
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APPENDIX G 


ADDITIONAL RESULTS FROM CONSERVATION EQUATIONS 


In addition to the conservation equations themselves, some interesting 
and useful side results can he obtained from the above derivation of the con- 
servation equations. The energy equation presented above, in its complete 
form, can be used to calculate the heat transfer and the temperature distribu- 
tion in the body, including the depth of the layer of molten stone and the 
depth at which the solid stone receives significant heating. 

For this purpose, equation (B26) may be written (neglecting P^/pg in 
comparison to unity) as 


" (lF) a = % + <3r + <Vb (}ab " RT b + 2 V) + (Pfa " T bK 

in which point a is arbitrary, so that 

^ab = c sol^f ~ ^a^ + ^f + c liq(-^b ^f) + ^v 9 ( x a — x f) 

^ab = c li<i( T b “ T a) + L v * < x f < x a < °) 

Then equation (Cl) becomes 

(dT/dx) - p 1 (T - Kj = 0 , (x < x f ) 

(dT/dx) - P 2 (T - K 2 ) = 0 , (x f < x < 0) 

where 

= p b u b c sol = Pb u b C liq 

' *sol ' ^ - kl . q 


(Cl) 


(C2) 


(C3a) 

(C3b) 


(C4) 


^i K i k sol \^b + ^-r + ^ u b c sol T f + L f + c liq( T b “ T f) 


+ L v - RT b + 2 U ^ 2 


+ (Pb “ T b) u b} 


(C5a) 


P2 K 2 = { ac b + Ir + Pb u b 


c liq T b + L v “ R^b + g ti^ 2 


+ (Pb - T b) u b f 


mill 


Much simpler expressions for K x and K 2 than equations (C5) can he 
obtained as follows: Take point a far enough back from the surface so 

q a is small and T a is the cold interior temperature. Then, as x -» -< 
T ~ T a and dT/dx 0; K x is found (from eq. (C3a)) to be 

K x = T a 

Further, equations (C4) through (c6a) can be combined to give 


K 2 = 


c sol Pi Tjr 

1 ' 

c sol 


Pb^b 


p 2 k 


liq 


C( G sol G liq)^f + ^f] 


= T f - 


: iiq 


(Tf - T a ) - 


Lf 

: liq 


The boundary condition for equation (C3*b) is 

x = 0 , T = T b 

Since T^ is known, the location Xf will be determined by solving 
equation (C3h) for the condition 

x = Xf , T = T f 

The boundary condition for equation (C3a) is then 

x = x f , T = Tf 

The solutions to equations (C3a) and (C3b) are 


T = Kx + K 1 'e P i X , 

(x < x f ) 

T = K 2 + K 2 'e P2X , 

(x|> < x < 0) 


Substitution of the boundary conditions gives 


T b = K 2 + K 2 t 

T f = K 2 + K 2 'e® 2Xf > 


Tf 


Ki + Ki'e 


PlXf 


J 


from which the expressions for K 2 r , K x r , and Xf are obtained: 


k-2 


■ ■■ an i i n n 1 1 in ii 


I 


that 

(c6a) 

(C6b) 

(CT) 

(C8a) 

(C8h) 

(C9a) 

(C9b) 

(CIO) 


I 


Note also that 


Ka' = T b - Ka 

_ _1_ . A f ~ Kg\ 

x f p 2 lo S _ k 2 J 


Kx* 


(Tf 



K 2 \-Pi/P 2 

K z) 


J 


Ko.* = (T f - Kje' 31 ^ 1 
K 2 ' = (T f - K 2 )e" p2Xf j 


(Oil) 


(C12) 


and that the distribution (C9) can be written conveniently as: 


where 


T ~ Ki Pi(x-xf) 
Tf - K X 

T ~ K 2 _ p 2 (x-Xf ) 
Tf - K 2 6 


(X < Xf ) 

> 

(xf < x < 0) 


(C13) 


Kx = 
K 2 = 

x f = 


m ISO 

T f - c Li 

it 108 ( 




_Lf_ I 
c liq ( 


J 


(ciiO 


and where Pj_ and p 2 are given by equations (c4) . Thus, the temperature 
distribution in the body is determined (eqs. (C13)), as well as the depth of 
the molten layer, -Xf (eqs. (Cl4) ) . The depth of the solid stone material 
receiving significant heating can be observed from the temperature distribu- 
tion, or one can calculate the depth at which the temperature is, say, 1 per- 
cent higher than the cold interior temperature. If one calls this depth -x c 
and the temperature T c , the first of equations (C13) gives 


Tc ~ Ki Pi(x c -x f ) 
Tf - Kx 


from which 

x c = Xf + p7 log (tJ - k£) 

where, for example, T c = 1.01 T a . 


& 
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To compute the temperature distributions, one needs to know values for 
k so l and kliq> in addition to the constants listed in table II. Opik gives a 
value for thermal conductivity of meteoric stone as 2X10 5 erg/cm sec °K 
(ref. 12, p. l62) • We take this to be 

k so l = 2xl0 5 erg/cm sec °K 

In general, the ratio kqj_q/k so ]_ is about one -half for a material at fusion. 
Therefore, we assume 

kfiq = IX10 5 erg/cm sec °K 

We find: K ± = 200° K, K 2 = 257° K, Pi/p^ = 44.75 cm sec/gm, and 
Ps/Pi = 2.46. For the Pfibram meteor, the depth of the molten layer and the 
depth for 1 percent temperature rise are shown in figure 8. The interior 
temperature distribution at 55 km altitude is shown in figure 9* The calcula- 
tions show that the depths of significant heating in the meteoric body are 
very small. 



Figure 8.- Depth of molten layer (-x f ) and depth 
of 1 -percent temperature rise from cold state 
(-x c ) (x-b* = 20 cm, s = 2). 



Figure Representative interior temperature 
distribution (55 km altitude, r^* = 20 cm, 
s = 2) . 



APPENDIX D 


CALCULATION PROCEDURE 


The flight parameters are calculated as follows: 

1. Specify as input: 

Tables of altitude and flight data: g sin 0 and h vs. p 

(table I); h vs. V and dV/dt (table V(c) for r^* = 20 cm) 

Constants: p gL (eq. ( 2 )); Mb (eq. (38)); s (eq. ( 21 d)); p m , c sol , 

c liq> T a , T f , L f , L v , R, and 7 (table II) 

2 . Specify altitude h 

3. Compute shock density ratio k = P^/Ps (eq. (3)) 

4 . Compute D (function of radius and standoff distances; see eq. ( 4 )) 
by using trial values for D and iterating until A(k, D) = 0 (eq. (6)) 


5 . Compute : 


P = Pot ■ P = ambient density at altitude h 
” SL (e 3 . («) 

p Q = p^V 2 ~ stagnation pressure (eq. (15)) 

I = l(p, V) = outer shock-layer radiation rate 

per unit volume behind normal shock (eq. (7)) 



(eq. (10)) 

^OO 

8 P m( -dV/dt + g sin 0) 

1 (Tv^ 


(eq. (13)) 


Z 



(eq. (IO)) 


6. Determine Mg and other quantities by using trial values of M e and 
iterating until a = Z - Z l = 0 ^ where the following quantities must be 
computed for each trial value of Kq: 
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P e = P e (p Q , M e^ 7) (eg., (l 6)) 
a = a(Mg, 1 ^, 7 ) (eg. (40)) 

P e /Pb = 1/a (eg. (39)) 

T e /T b = function of 7 , a (eg. (4l)) 

P e /Pb = aM b 2 /M e 2 (eg. (42)) 

Ch c = function of a, 7, M e , ^ (eg. ( 43 )) 

_ p e 
^ Pe/Pb 

T b = T t>(Pb) (e<1 - ( 36b)) 

^ab = c sol(T f - T a ) + L f + c li g(l'b “ T f) + p v ( e< l- (25)) 


Z' = Z»( 7 , R, T b , Mb, £ ab , C hc ) (eg. (35)) 
cj = Z - Z* 

7 * Then compute: 

Pb = Pb/RTb ( e< l- (34)) 
u b = M b N^ib 


(chJ 


r f f 
c H r 


2pb u b 3 Aab A , 

grv^A + 2) (e *- P° 


s PH r ) ff 


r b " 


(CH r ) ff 


(CH r /^b) 


ff 


C D. 


eff y 


D eff 


Ph u b 

(CH c ) ff = C hc (eg. (31)) 
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CH c = I ( c H c ) ff (ei- (50)) 


Ct) = - 
•^p s 


Pb 


(eq.. (47)) 


CD- 


2 Poo V " 

V(CHc)« 

ff - (eq. (U8)) 


VC 

su b 

Pt = 

°Peff 


b7 


APPENDIX E 


ESTIMATION OF THE VAPOR -ABLATION REYNOLDS NUMBER 


The validity of the equations for the boundary shock wave (ref. l) 
depends on the vapor -ablation Reynolds number 


Re b = 


_ frbVi 


P b 


(El) 


being large (» 1) . The quantities p^ and u^ are evaluated from calcula- 
tions presented in the text. To obtain a rough estimation of Re^, in order 
to determine rough limits on the applicability of the theory, one must 
estimate dj_ and p^ . 


The shear viscosity, q, to which p is related by 


P = 



(E2) 


can be estimated using approximate equations from kinetic theory (e.g., see 
Vincenti and Kruger, ref. 7 )• Let w be the molecular weight of the vapor 
and consider the vapor to be composed of fictitious identical molecules of 
diameter d. Let be the mass density per unit volume of the liquid state 

of the vapor whose molecular weight is ^ . Then, denoting Avogadro T s number 
by £, where 


£ = 6 . 023X10 23 


molecules 

mole 


we may write for the approximate volume occupied by each molecule in the 
liquid state 


3 _ 2 % 


d = 


p z i 


(E3) 


The coefficient of shear viscosity (see ref. 7) may be approximated by 


1 K c 

^ 2 * J2 £ a 2 


(E4) 


where 


C« (E5) 

Substitution of equations (E 3 ) and (E5) into (E4) gives the approximate 
relation 
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a/ 3 


n - ^GsLAVihf 

2 J2 rt \ZJ\M 


Now If one assumes 


P 


l 


P 


m 


and takes values for R, and p from table II. one obtains 

m ? 


-JI 

1 gm cm*" 1 sec"* 1 


o. 253x10 “ 4 



(E6) 


(E7) 


(E8) 


For convenience, one may assume the bulk viscosity k to be zero, so that 
(from eq. (E2) ) 

P = 3 4 (E9) 

Then, with the result from equation (38) that 

u b = O.333 N//RTb (E10) 

equation (El) with (E8) and (E9) becomes 


Re b = 


fib u b d i 

h 


10 7 ( fib \ / r b \ (*£\ 
0.498 VI g ra cm“ 3 y \1 cm J V r b/ 


(Ell) 


The estimation of d b /r b lies outside the scope of this report. It has 
been computed from an approximate relationship for which the derivation is 
rather involved. (The flow field around a sphere with large mass transfer was 
analyzed.) It will simply be noted here that, corresponding to the calcula- 
tions for the Pribram meteor (with r^* = 20 cm and s = 2) , the following 
approximate values for di/r^, and the resulting Reynolds numbers, Re b , were 
found : 


h 

a i 

Re b 

1 km 

r b 

75 

0.00507 

0 . 204 X 10 

70 

.00939 

.I38XIO 2 

65 

.01629 

.767XIO 2 

60 

.02661 

. 354x10 s 

55 

.03928 

. 120X10 4 

50 

.05039 

• 268 X 10 4 

45 

.03407 

.857x10 s 


h-9 



With little extra effort one can also estimate the thermal conductivity 
in the vapor: 


V ^°P 4 4 (L T) 

k = ;^r~ « ~ [ic tv 7 \i ( ir R 
Pr 3-^3 


14 


I \±R 


(E12) 


With equation (E8) and the value of R from table II, equation (E12) gives 

>1/2 


, Qfi / ._T 

1 gm cm sec _ 3 (°K) -1 y l 1° K 


(E 13 ) 
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TABLE I . - TRAJECTORY DATA FOR ONDREJOV METEOR PRIBRAM AMD 
ALTITUDE -DENSITY FUNCTION FROM REFERENCE 3 
sin 0 = O.6853; f = g sin 0 = 671.5 cm/sec 2 


Altitude, 

Density function, 

h/1 km 

P 

90 

0.284, -5 = 0. 28 4x10 ~ 5 

85 

.634, -5 

80 

• 173, -4 

75 

•393. -4 

70 

.817, -4 

65 

.158, -3 

60 

.288, -3 

55 

•499, -3 

50 

.884, -3 

45 

.164, -2 


TABLE II . - APPROXIMATE VALUES OF PHYSICAL CONSTANTS 


APPROPRIATE TO A STONE METEOR 


Symbol 

Property 

Value 

Unit 

Reference 

Pm 

Mass density of meteor 

3-5 

gm/cm 3 

3 

1 — 1 

O 

W 

O 

Specific heat in solid 

8.95x10 s 

erg/gm °K 

12 

c liq 

Specific heat in liquid 

1. ixio 7 

erg/gm °K 

12 

T a 

Typical cold interior temperature 

200 

°K 

12 

T f 

Fusion temperature 

1800 

°K 

12 

Lf 

Latent heat of fusion 

2.65x10 s 

erg/gm 

12 


Latent heat of vaporization 

6.05X10 10 

erg/gm 

12 

-r?z 

Molecular weight of vapor 

50 

gm/mole 

_ ’ | 

12 

R 

Gas constant 

l. 66xio 6 

cm 2 /sec 2 °K ! 

3 

7 

Ratio of specific heats in vapor, ' 

C p/ C V 

1.4 


1 
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TABLE III.- METEOR -TRACKING DATA FOR ONDREJOV 
METEOR PRIBRAM (Ref. 4 , table 7) 


t, sec 

h/l km 

V/lO 5 cm sec " 1 

0 

.85806 s 

.85602 s1 

1.73230 

2.49383 

2.69203 

3.06758 

88 . 594 

76.289 s 

76.318 s 

63.837 

52.970 

50.164 

44.858 

20.887 ±0.009 
20.864 ±.oio a 
20.860 ±.007 a 
20.838 ±.013 
20.773 ±-013 
20.717 ±.013 
20.459 ±.024 


& The data at these two points 
(both at about h = 76 km) were from 
two different rotating -shutter 
cameras. The first two points on the 
table are from one camera, the 
remaining five points from the other 
camera. 


TABLE IV.- CALCULATIONS IN CONSTRUCTION OF LATA CURVE 


h 

1 km 

Values of 

\10 5 sec 

^ for r-^" 

* equal 

to 


10 cm 

15 cm 

20 cm 

25 cm 

35 cm 

50 cm 

90 

-35.6 

-39.4 

-41.3 

- 42.45 

-43.65 

-44.7 

85 

-21.65 

-30.08 

-34.35 

-36.8 

-39.7 

-41.9 

BO 

22.03 

-.958 

-12.4 

- 19.35 

-27.2 

-33.1 

75 

109.7 

57.5 

31.5 

15.8 

-2.16 

-15.5 

70 

279.0 

170.3 

116.2 . 

74.0 

46.15 

18.45 

65 

584 

373 

269 

205.3 

133.2 

79.5 



TABLE V. - VELOCITY AND ACCELERATION DATA FROM GRAPHICALLY FITTED CURVES 


h 

V 

dV/dh 

-dV/dt 

1 km 

10 5 cm sec -1 

10“ s sec -1 

10 3 cm sec -2 


(a) r t 

* = 10 cm 


90 

20.877 

-35.6 

-0.510 

85 

20.878 

-21.65 

-.310 

8o 

20.8775 

22.03 

.3155 

75 

20 . 874 . 

109.7 

1.569 

70 

20.863 

279.0 

3.985 

65 

20.842 

584 

8.350 

6o 

20.811 

743 

10.60 

55 

20.767 

1,085 

15.43 

50 

20.688 

2,255 

31.97 

45 

20.468 

8,330 

116.95 


(t>) r^* = 15 cm 


90 

20.860 

- 39.4 

-0.564 

85 

20.861 

-30.08 

-.430 

80 

20.862 

- 0.958 

-.0137 

75 

20.860 

57-5 

.8225 

70 

20.853 

170 

2.430 

65 

20.841 

373 

5.340 

60 

20.820 

527.5 

7.530 

55 

20.783 

960 

13.670 

50 

20.710 

2, 160 

30.64 

. 45 

20 . 470 

9,800 

137.4 


(c) 17 

* = 20 cm 


90 ' ' 

20.859" 

-41.3 

-0.590 

85 

20.861 

-34.35 

-.4905 

80 

20.862 

-12.4 

-.1773 

75 

20.862 

31.5 

.450 

70 

20.8575 

116.2 

1.66 

65 

20.848 

269 

3-84 

60 

20.8285 

533 

7.61 

"55 

20.791 

1,023 

14.56 

50 

20.714 

2,300 

32.63 

45 

20.485 

10,400 

145-8 


(d) 1^* = 25 cm 


90 

20.8585 

- 42.45 

-0 . 606 

85 

20.861 

-36.8 

-.526 

80 

20.862 

- 19.35 

-.2765 

75 

20.8615 

15.8 

.226 

70 

20.860 

74.0 

1.059 

65 

20.853 

205.3 

2.933 

60 

20.8375 

445 

6.360 

55 

20.802 

1,021 

14.550 

50 

20.711 

2,910 

44.100 

45 

20.465 

7,150 

100.2 


55 



TABLE V. - VELOCITY AMD ACCELERATION DATA FROM GRAPHICALLY 


FITTED CURVES - Concluded 


h 

V 

dV/dh 

-dV/dt 

1 km 

10 5 cm sec~i 

10 -5 sec -1 

10 3 cm sec"* 2 


(e) ] 

r b * = 35 cm 


90 

20.854 

-43.65 

-0.624 

85 

20.8565 

-39.7 

-.567 

80 

20.858 

~27.2 

-.3885 

75 

20.8585 

-2.16 

-.03085 

70 

20.8575 

46.15 

.660 

65 

20.852 

133.2 

1.903 

60 

20.837 

48o 

6.855 

55 

20.800 

1080 

15.4 

50 

20.711 

2750 

39.0 

45 

20.4795 

8340 

117.0 

(f) rt* = 50 cm 

90 

20.852 

-44.7 

-0.639 

85 

20.854 

-41.9 

-.599 

80 

20.856 

-33.1 

-.473 

75 

20.8575 

-15.5 

-•2215 

70 

20.857 

18.45 

.264 

65 

20.854 

79-5 

1.137 

60 

20.843 

373 

5.325 

55 

20.812 

934 

13.300 

50 

20.724 

2680 

38.00 

45 

20.480 

8070 

113.4 


56 
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